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Abstract 

Recently the NA60 collaboration has reported the transverse mass spectra of dimuons coming 
from In- In collisions at 158 GeV/A. The measured yields display a strong invariant mass de- 
pendence not typical of radial flow, suggesting that different sources contribute in different mass 
regions. We interpret the dimuon transverse mass spectra from an early thermalized partonic 
phase and hadronic phase constrained by the strictures of broken chiral symmetry. Each phase 
develops a specific transverse momentum dependence by hydrodynamical expansion. We show that 
a measurement of the momentum anisotropy at NA60 could provide information on the dominant 
emission source (hadronic or partonic) in the intermediate mass region 1.5 < M < 3.0 GeV. 

PACS numbers: 
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I. INTRODUCTION 



The observed dileptons coming from heavy-ion collisions can provide useful information 
on the early stages of the collision since electromagnetic probes do not interact with the 
medium after they are produced. It is therefore hoped that observed dileptons and photons 
could possibly provide an unambiguous signal for the quark gluon plasma. However, in 
any collision, there is also a substantial contribution of dileptons coming from the hadronic 
phase. An understanding of the resulting hadronic yields can provide crucial information 
on modifications to electromagnetic spectral functions due to chiral symmetry restoration. 

The recent NA60 experiment at the CERN SPS has measured the invariant mass spec- 
trum of low-mass dimuon pairs in In-In collisions. It was seen that a large excess 
remained after subtracting contributions from expected hadronic [the cocktail) decays. The 
remaining excess was examined by a number of groups 0, 0, Q, s[ [9] and was interpreted 
by a combination of thermal partonic and hadronic contributions with modifications to the 
spectral function due to finite temperature and baryon density. 

Even though the spectral function directly gives the rate of dilepton production, the 
experimentally observed yield differs due to the fact that the rates must be convoluted 
over the full space-time history of the collision region having widely varying temperature, 
baryon density, chemical potential and flow velocity. In addition, contributions from an 
early partonic phase must be included. 

The NA60 results for the invariant mass spectra are important for pinning down modifi- 
cations to the spectral function at finite temperature and baryon density but do not provide 
any information on collective transverse flow as the spectra are Lorentz invariant (not taking 
into account the nontrivial acceptance at NA60). More recently the NA60 collaboration has 
released acceptance corrected transverse mass spectrum 4| of the resulting dimuon yields in 
different mass windows. It was found that the spectrum has a shape atypical of radial flow 
with a significant mass dependence suggesting that the different mass regions are populated 
by different sources. 

In the following paper we continue the analysis in {5] and analyze the transverse mass 
spectra of the dimuons produced at NA60 using the same rate equations and hydrodynamic 
evolution. In section 2 we summarize the pertinent partonic and hadronic rates used. In 
section 3, we present our transverse spectra and elliptic flow predictions. Our discussion 
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and conclusion follows in section 4. 



II. DILEPTON RATES 

We now summarize the rate equations used for the analysis that was preformed in |5]. For 
the partonic contribution above a critical temperature T c « 170 MeV we use the standard 
leading order qq result for massless quarks: 



dR 
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12vr 4 ei' 



,71-1 
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(1) 



q=u,d,s 

where Nc is the number of colors, e q the charge of t he q uarks, and n± = l/(e^ 0± '^^ 2r + 1). 



Below T c we use the rate equations presented in 



to 
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1 2t ] for a hadronic gas at finite 



temperature and baryon density which are constrained entirely by broken chiral symmetry 



dR = -a 2 1 2m 2 _ 4mf 1/2 

d 4 g 3vr 3 g 2 1 + e« / T 1 g 2 Jl g 2 J 

x [-3g 2 ImII y (g 2 ) + ^ / ^Wf (g, fc) + ^ ciiVW^(g, p)] (2) 

where da and dN are the appropriate phase space factors for mesons and nucleons respec- 
tively as outlined in [5] . The term in square brackets is obtained after keeping terms to first 
order in an expansion of the thermal structure function in both meson and nucleon density. 
Through the use of three-flavor chiral reduction formulas the terms Wf (g, k) and Wjy(g,p) 
can be expressed in terms of vacuum vector and axial spectral densities which are measured 
in e + e~ annihilation, r decays and photo-reactions on nucleons and nuclei. 

The final contribution comes from the decay of vector mesons after freezeout. These 
dileptons will play a large role in the transverse momentum spectra in the mass region of 
the given meson. The main effect is to harden the transverse mass spectra as the decaying 
mesons at freezeout are coupled to the transverse flow of the medium and are therefore 
boosted to higher q?. Since the contribution from the decay of freezeout vector mesons was 
not included in our first analysis Q], we discuss the rate equation and the integration over 
the freezeout surface in the appendix. 
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FIG. 1: Resulting dimuon yields after integration over the space time volume of the collision 
convoluted with a schematic acceptance of the NA60 experiment. Shown are the contributions 
from the perturbative QGP and charm decays (pQGP + DD), the decay of freezeout p mesons 
(FO), and the thermal hadronic yields (Hadronic). 

III. RESULTS 



The rates outlined in the previous section were integrated over the space-time history of 



a full hydrodynamic simulation (details discussed in [5|) of the collision region at the CERN 
SPS collider. Shown in fig. [I] are the resulting invariant mass spectra (after a schematic 
acceptance of the NA60 detector was applied) compared to the data measured by the NA60 
collaboration. The only difference from the analysis done in |5[ is the inclusion of the 
freezeout contribution as shown in the figure. The net result is the same within %10 except 
for a small change in the overall normalization of the rates which can be accounted for by a 
decrease in the fireball freezeout temperature. 

The data points in fig. [2] show the acceptance corrected transverse momentum spectrum 
of observed dimuon pairs after subtraction of the known hadronic cocktail (omitting the p) 
and charm contribution. The data was divided into three mass windows: a low mass region 
(0.4 < M (GeV) < 0.6) below the p peak, a mass region around the p (0.6 < M (GeV) < 
0.9) and also a higher mass region (1.0 < M (GeV) < 1.4). The data selection contains 
dN c h/dr] > 30 which corresponds to about < dN c h/dr] >~ 140 which we use in the analysis. 

The solid lines in fig. [2] show the results using the rates discussed in the prior section. 
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We show the contribution to the qr spectrum from three sources: partonic (QGP), hadronic 
(Had) and freeze-out (FO) p contribution. The upper line in each figure shows the total 
yield. 

Since a full hydrodynamic simulation of the collision is available we go one step further 
and calculate the elliptic flow of the produced dileptons [16( as is shown in fig [3] where v2 
is a function of the invariant mass, transverse momentum and the impact parameter of the 
colliding system. The elliptic flow is defined as the weighted average of the yields with 
cos(20): 

V2 ^ qT) = J C ° S(20) dMMyq T dq T d^ '/ W dM>dZ*tr*l> (3) 
where <fi is the angle between the dilepton momentum (q) and the corresponding fluid ele- 
ment's velocity. 



IV. DISCUSSION AND CONCLUSION 



Our rates are overall consistent with the data in all mass regions. There is the possibility 
of some enhancement in the high q? regions in both the low mass and medium mass regions. 
This could possibly be accounted for by fine tuning the freezeout temperature thereby in- 
creasing the freezeout p contribution and thus hardening the spectra in these mass regions. 
The drawback however, is that our yield around the rho is already about 30% above the 
measured yield. There is also some discrepancy in the low qx regions, but as was pointed 



out in 



201 ] . one should not expect to explain this region by a thermal description only. 



Based on the data, we feel that overall we have good control over the hydrodynamic 
evolution of both the temperature and flow velocity as well as the different source contribu- 
tions (hadronic or partonic) in different mass regions. Of course with additional fine-tuning 
of fireball parameters (freeze-out temperature, etc.) it is possible to achieve better agree- 
ment, but we don't feel this is worth doing since the data is integrated over a wide range of 
centralities. 

There has also been discussion whether the mass region above 1 GeV is dominated by 
partonic or hadronic radiation. Even though information on the source can be extracted 
from q? spectra by looking at the effective temperature of the emission region as was done 
in [2^] more information could be obtained if it was possible to measure the dilepton elliptic 
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FIG. 2: (Color Online) Transverse momentum spectra for the three mass bins as measured by 
NA60. From top to bottom: 0.4 < M (GeV) < 0.6, 0.6 < M (GeV) < 0.9, 1.0 < M (GeV) < 1.4. 
The solid lines are our results including contributions from the partonic phase (QGP), hadronic 
phase (Had) and vacuum decays after freezeout (FO). 
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FIG. 3: (Color Online) Dilepton elliptic flow as a function of qx for semi-central In-In collisions. 

flow in various mass windows. Since the thermal partonic radiation is emitted early on from 
the collision region, one does not expect a large momentum anisotropy to have developed 
and therefore one would expect a much lower elliptic flow from dileptons coming from this 
region. On the other hand, hadronic radiation is emitted much later on and any momentum 
anisotropy would also be carried by the emitted dileptons. 

In our rate calculations we find that in the 1.0 < M (GeV) < 1.5 region the dilepton 
contribution comes from both the hadronic and partonic phases. The hadronic yields in 
this mass region comes from contributions of both the axial spectral density which enters 
the rates from the chiral reduction formula due to chiral symmetry restoration by matter 
and also from the high energy tail of the vector spectral weight, part of which is from 4-7T 
annihilation. In fig. [3] we show our calculated v2 from the partonic and hadronic yields. 
We also show the total v2, which is a weighted average of the two sources. One can see 
that depending upon the dominant emission source that qualitatively different elliptic flow 
develops at larger qx- 

To summarize: We have interpreted the dimuon transverse mass spectra as measured 
by the NA60 collaboration as coming from an early thermalized partonic phase as well as 
a hadronic phase constrained entirely by broken chiral symmetry, both of which develop 
specific transverse momentum dependence following the underlying hydrodynamic flow. 
After convolution of the rates over a hydrodynamic simulation of the conditions at the 
CERN SPS experiment we find good agreement between the resulting yields and the 
measured data. We show that a measurement of the momentum anisotropy at NA60 



7 



could provide information on the dominant emission source (hadronic or partonic) in the 
intermediate mass region 1.5 < M < 3.0 GeV. 
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APPENDIX A: DILEPTONS AFTER FREEZE-OUT 



We discuss how the dilepton rates from vector mesons decaying after freezeout was evalu- 
ated in the above result. A number of works [l7, isl, 19] discuss this topic but it is presented 
here as well for completeness. 

The contribution to the dilepton spectrum from decaying vector mesons at freezeout is 
given by: 

dN _ y 1 dTy^^^- d 3 N v A _ 
dM 2 d 3 q/q° ^ Vf dM 2 d 3 p [P q> { ' 

In the above equation p is the momentum of the decaying vector meson V, q = fa + fa is 
the momentum of the dilepton pair and M = + is the invariant mass of the dilepton 
pair. The quantity E^- is the number of vector mesons with momentum p and energy E 
at freezeout and can be evaluated using the Cooper-Frye [iJ] formula 

E ^ = 7^ f f(qo,T)p»da, (A2) 



d 3 p (2tt) 3 

The Cooper-Frye formula requires integration over a, which is the four-dimensional space- 
time boundary. We show how this integration was done for a boost invariant, azimuthally 
symmetric expansion. 

The quantity da^ is the normal vector to the surface of the four- dimensional space-time 
volume set in our case by the freeze-out temperature 7/. . and is given by 

, da v da x daP , , , . , 

d a M = tfiuXp q d^j- du dv dw (A3) 

where u,v, w are three independent coordinates used to parameterize the hypersurface. Fol- 



lowing 



151 ] the four volume of the collision region is expressed as 



<7 m (t, r, (p, rj) = [t cosh rj, r cos0, r sin0, r sinhr/] (A4) 



where r\ = | In |^ is the spatial rapidity and r = \Jt 2 — z 1 is the proper time. Also for a 
boost invariant expansion, the vector meson's four momentum can be expressed as: 

]f = [m± cosh y, p±_ cos <fi p , p± sin P , sinh y] (A5) 



where m_i_ = a/ M 2 + an d y is the longitudinal rapidity. 

Choosing to parameterize the freezeout hypersurface with coordinates (u,v,w) = (r,0,?7) 
the freezeout volume is: 



Ot\r. 



V) = [ T f ( r ) cosh 77, r cos 0, r sin 0, Ty (r) sinh 77] (A6) 



The subscript f in the above equation refers to the quantity at freezeout. Therefore 
should be thought of as a parameterization of the space-time volume where freezeout occurs, 
or more precisely in our case the space-time volume having temperature T = Tf_ Q _. The fact 
that the freezeout proper time depends only on the fireball radius is due to the assumed 
boost-invariance and azimuthal symmetry. Using lA3l we obtain 

d 3 Oa = (cosh ti, — — — cos 0, — — sin 0, — sinh rifTrdrdcftdri (A7) 
^ Or or ' 

for the normal vector to the hypersurface. 

Using the above expressions for p^ and da^ and performing the integration of Eq. IA2l over 
the freezeout hypersurface, the total number of dileptons produced from decaying p mesons 
per unit four momentum is given as: 

dN n°- , f +QO , f 2w 71 
— — = dr drj I dcp rr 

d Q Jo J-00 Jo 

x [m ± cosh( 27 - rj) - %p ± cos(0 - 0^)] ° 2 T) (A8 



dr ^ T JS Tf dM 2 (2 



7T 



where we have chosen to orient the azimuthal angle with the direction of the fluid 
velocity and f{p ll u ti ) T) is the Bose distribution with p^u^ = m±'y± cosh(y— rj) —p±uj_ cos(0— 
<ff). The explicit form for the meson spectral function is given using the Breit-Wigner 
parameterization: 
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dT 



dM 2 



7T (M 2 - ml) 2 + (mpTfty 



(A9) 
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